Evolution of the correlation functions in 2D dislocation systems 
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In this paper spatial correlations of parallel edge dislocations are studied. After closing a hierarchy 
of equations for the many-particle density functions by the Kirkwood superposition approximation, 
we derive evolution equations for the correlation functions. It is found that these resulting equations 
and those governing the evolution of density fields of total as well as geometrically necessary dislo- 
cations around a single edge dislocation are formally the same. The second case corresponds to the 
already described phenomenon of Debye screening of an individual dislocation. This equivalence of 
the correlation functions and screened densities is demonstrated also by discrete dislocation dynam- 
ics simulation results, which confirm the physical correctness of the applied Kirkwood superposition 
approximation. Relation of this finding and the linear response theory in thermal systems is also 
discussed. 
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I. INTRODUCTION 
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It is well known that plastic deformation of crystalline 
materials is caused by the motion of a large number of 
interacting dislocations. Although the properties of in- 
dividual dislocations have been known for decades, the 
question how to build up a micron-scale level continuum 
description of dislocations is still an unresolved issue in 
the theory of crystal plasticity. The appropriate contin- 
uum theory could provide a framework for understanding 
dislocation pattern formation or size-effects. 

Several phenomenological models have been proposed 
for dislocation patterning so far. Here we first men- 
tion the work of Walgraef and Aifantis^'^, who adopted 
a reaction-diffusion model from the description of oscil- 
lating chemical reactions. In their model second-order 
gradient terms are introduced to account for the spa- 
tial fluctuations of the dislocation density. Although the 
framework turned out to be very successful in the pre- 
diction of different dislocation structures, the physical 
origin of the length scales introduced through the coef- 
ficients of the gradient terms is not really understood. 
The same problem applies to the model of Kratochvil 
and co-workers, in which nonlocal terms are introduced 
for describing spatial interactions through the sweeping 
mechanism i^i^ 

In a non-phenomenological model proposed by Groma 
a system of parallel edge dislocations in single slip was 
studiedi^ By performing ensemble averaging on statisti- 
cally equivalent dislocation distributions a chain of equa- 
tions for the many-body dislocation density fields was 
derived. First, dislocation correlations were neglected, 
which corresponds to a mean- field approximation.^ Here 
the stress acting on dislocations is simply the "self- 
consistent" field, which is the long-range stress field of 
the geometrically necessary dislocations. If the disloca- 



tion system is correlated, which is indeed the case for 
real systems (see, e.g., Ref. 0), new stress-like terms ap- 
pear in the constitutive equationsi^ The most important 
non-trivial term is the gradient-like "back-stress" , which 
can be interpreted as the short-range effect of disloca- 
tion pile-ups. In order to check the validity of this two- 
dimensional (2D) theory, its predictions were compared 
with discrete dislocation dynamics simulation results and 
good agreement was foundi^iiiS 

There are several current approaches in the literature 
to develop a density-based continuum model for three 
dimensions (3D). El-Azab based his theory on Nye's dis- 
location density tensor. After studying the mean-field 
approximation of an uncorrelated ensemble, in their re- 
cent paper El-Azab et al. turned to the investigation 
of dislocation correlations on simulated 3D dislocation 
ensembles in a body-centered-cubic crystalJ^ii^ Csikor 
et al. reported about a complementary numerical study 
focused on the range of 3D dislocation pair correla- 
tions in face-centered-cubic materialsJ^ Recently, an- 
other promising theory has been proposed by Hochrainer 
and co-workers who studied the evolution of dislocation 
lines in a higher-dimensional space-i^iil Motivated by 
the 2D results dislocation correlations are introduced 
into the theory semi-phenomenologically by considering 
short-range stress terms such as the aforementioned back- 
stressJi Similar extension was adopted by Schwarz et 
al. tooi^ for a continuum theory of interacting curved 
dislocations^?. In all mentioned current 3D continuum 
theories of dislocation dynamics one of the biggest chal- 
lenges is the precise incorporation of the effects of dislo- 
cation correlations, which seems to be unavoidable. So, 
the investigation of the evolution of the correlation prop- 
erties is of great importance, but even in 2D this problem 
has not been solved yet. 

In the past two years new light has been shed on equi- 
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librated dislocation systems by the discovery of the ef- 
fect of dislocation screeningi22i21 Its most direct physi- 
cal significance is that it explains the extensivity of the 
elastic energy in equilibrium by the appropriate relative 
arrangement of the dislocations, resulting generically in 
a finite interaction length rather than the unscreened in- 
finite range interaction. On the technical side, in the 
simplest example of single slip, 2D, edge dislocation sys- 
tems, an effective thermodynamic potential was proposed 
for the purpose of the variational calculation of the ge- 
ometrically necessary density. Summarizing the results, 
for large and constant total dislocation density and small 
but variable geometrically necessary density, the equa- 
tion for the induced density with arbitrary external dis- 
locations was constructed and its Green function for in- 
finite boundary conditions was given analytically. Hence 
the induced density and the resulting screened elastic 
potential by a single external dislocation was obtained, 
showing suppression in all directions exponentially strong 
except for one axis, where it was of power type. This phe- 
nomenon bears close resemblance to Debye screening of 
an external charge in Coulomb systems, so this term was 
adopted also for dislocations. In the same work^° numer- 
ical results were presented, but because of the prohibitive 
fluctuations of the dislocation distribution induced by an 
external, fixed dislocation, rather the signed pair correla- 
tion of dislocations in equilibrium was monitored. Even 
in the latter function an agreement with the theory was 
found. 

It should be borne in mind that the relation between 
the externally induced density (response function) and 
the correlation function is based on linear response the- 
ory in thermal systems. However, in the dislocation sys- 
tem considered presently entropic effects are suppressed, 
i.e., thermal fluctuations are negligible beside the Peach- 
Koehler forces. Therefore, traditional arguments of lin- 
ear response theory do not obviously apply. This raises 
the question why then the externally induced dislocations 
obey properties similar to those of correlation functions. 
The problem can be actually posed in a broad sense, 
namely, to what extent the response function to an ex- 
ternal particle and the two-particle correlation function, 
in or off equilibrium, are similar. The evolution equa- 
tions for the one-point densities, with some simplifying 
hypotheses, are known, ^ and they have been linked to the 
variational approach by a phase fleld constructions!. The 
evolution equations for the two-point correlation func- 
tions, however, could not be brought to any treatable 
form so far. 

As a recent development, Vinogradov and Willis pub- 
lished an evolution equation for the correlation function 
of long parallel screw dislocations with the same Burgers 
vectors <^2i2^ They even solved the resulting equation for 
the static case analytically. It was found, however, that 
the generalization of this result for the case of two possi- 
ble Burgers vectors (positive and negative dislocations), 
or for edge dislocations, leads to a solution not decaying 
at infinity, thus cannot be correct physically. 



In this paper we investigate the evolution of correlation 
functions of long parallel edge dislocations with single 
slip. We start our analysis with the Bogoliubov-Born- 
Grccn-Kirkwood-Yvon (BBGKY) hierarchy, describing 
the evolution of dislocation many-body densities, which 
was derived earlier by Groma. - In contrast to the clus- 
ter approximation used by Vinogradov and Willis ^^i^ 
we apply the Kirkwood superposition approximation as 
a method to close the chain of equations. The resulting 
equations, when rewritten for appropriately introduced 
single-point fields, are similar to those describing the evo- 
lution of one-dislocation density functions in the pres- 
ence of a single external dislocation. In fact they become 
identical, if the external perturbation is small enough, 
thus we arrive at the analog of linear response theory 
in thermal systems, now obtained for density fields at 
zero temperature. It should be kept in mind, however, 
that the relation is now conditioned upon the Kirkwood 
approximation for the two-point correlations. In order 
to test the extent of the analogy between the one- and 
two-point densities, we also performed extensive simula- 
tions. First we compared the equilibrium distributions 
induced by the screened single external dislocation with 
the appropriate fields from the two-point correlations, 
and found a quite satisfactory agreement. Second, the 
time evolutions of the two types of density fields were 
matched and again close similarity was observed even be- 
fore reaching equilibrium. This essentially demonstrates 
that the Kirkwood closure approximation was justified 
and the resulting equations indeed well describe the evo- 
lution of the correlation functions. We also analyze the 
cluster approximation of Vinogradov and Willis, adopted 
for edge dislocations, and show that it leads to a phys- 
ically unacceptable result for the evolution equation of 
the correlation functions. 

The paper is organized as follows. In Sec.|TT]we derive 
the new equations of motion for the two-point correlation 
functions and for the formally introduced, effectively one- 
point, fields. Section IIIII recuperates the equations for 
the single dislocation densities, while Sec. IIVI is devoted 
to the comparison of the previous two sections' results. 
Section |V] contains the presentation of the outcome of 
the simulations. The Conclusions are followed by the 
Appendix, containing background calculations for Sec. 

m 



II. TIME EVOLUTION OF THE 
CORRELATION FUNCTIONS 

Let us consider a system of parallel edge dislocations 
with single slip system, where only overdamped glide mo- 
tion is allowed. By considering a cross-section of the 
crystal perpendicular to the dislocations the problem be- 
comes two dimensional with point-like objects. If Si de- 
notes the sign of the ith dislocation, then its Burgers 
vector is bi — Sib [si G — } and b :— (b, 0)]. 

In order to describe a discrete system of dislocations, it 
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is useful to introduce discrete n-particle dislocation den- 
sities, e.g., jof ''*(r,t) := YliLi Ss,siS{r-ri{t)) forn = 1 in 
a system of N dislocations (for details see, e.g., Ref. I21I ). 
Here ri{t) denotes the position of the ith dislocation, and 
with the superscript s we make a distinction between the 
densities of positive and negative dislocations. Higher 
order discrete densities could be defined accordingly. In 
order to get smooth density functions, one has to perform 
averaging over statistically equivalent systems, which is 
denoted as, e.g., pl{r,t) := (pf''*(r, i)). After perform- 
ing the averaging procedure, Groma obtained a chain 
of evolution equations for the p*i' - '*" n-particle density 
functions of the system^S, These equations give the rate 
change of p^^' "'*", which depends on both itself and the 
one-level higher density function pf/^'x '^"^^ . This type of 
chain of equations is called BBGKY hierarchy in statis- 
tical physics. To make it possible to solve these infinite 
number of equations the hierarchy must be chopped at a 
certain level. 

We start our investigations at the first level of the hi- 
erarchy, which governs the evolution of the one-particle 
density functions!^ 

dt 



this case the system is homogeneous, meaning the one- 
particle density functions pi cannot have a spatial depen- 
dence, i.e.. 



:= p^{r,t) ~ const. 



(5) 



This is a trivial solution of Eq. ([T]). One can draw two 
simple conclusions from it: 

• The two-particle density functions depend only on 
the difference of their arguments: 

P2''^iri,r2,t) = P2''^iri -r2,t), si,S2 e {+,-}, 

(6) 



• and the following identity holds: 

^ S2 /" pl"''{r,t)nnd{r)d\ = 0, si G {+,-}. 

. Li 



S2-±l 



(7) 



We now proceed to the evolution equations of the two- 
particle density functions'!^ 



dpl'''^r^,r2,t) 
dt 

d d 



^''1 LI J Pi r 



[p^^^^^(ri,r2,t)F^-^Mn-r2) 



S2=±l 



+^[p?(n,i)Fex\(n)] = 0, 



(1) 



where p*^' ' '*" denotes the n-particle dislocation density 
function with Si referring to the sign of the ith disloca- 
tion, F'^^'^'^{ri — r2) is the interaction force between two 
dislocations at positions ri and r2, and with signs si and 

S2 (n eM^): 



^ E /'p3'''''''(n,r2,r3,i)F^-^nr2-r3)d2r3 



'■{r) = siS2Tindir)b. 



(2) 



+ Q^^[P2'^'^r,,r2,t)F:^,ir,)] 

+ ^[P2""(n,r2,t)F<f^,(r2)]^0. 



(8) 



Here Tjnd is the shear stress field of a single edge disloca- 
tion positioned in the origin: 



Tindir) = 



27r(l-iy) (a;2 + y2)2 



(3) 



in which /i and v denote the shear modulus and the Pois- 
son's ratio, respectively. A Text('') spatially varying ex- 
ternal shear stress applied to the system results in the 
extra F^^^^{r) force, which depends on the sign s of the 
dislocation: 



(4) 



Solving ([T]) generally (with arbitrary boundary condi- 
tions) requires the additional knowledge of p2 ■ 

In this section we investigate the relaxation of disloca- 
tions from a random initial state in an infinite medium 
at zero external shear stress. Since the equations of mo- 
tion of the dislocations don't have a direct spatial de- 
pendence, the system is invariant under translations. In 



It has to be noted that by definition P2^'^^ (ri, r2, = 

P2""(r-2,n,i). 

A possible method to truncate the hierarchy at this 
level is to express pg^'*^''*^ in the so-called cluster expan- 
sion, and neglect the three-particle correlation term, as it 
was done by GromE*^, and more recently by Vinogradov 
and Willi a^^'^'^ . This means pg^'*^'*^ is approximated as 

P3"'"'^(n,r2,r3,t) 

«pr(ri,i)pf (r2,t)pf (r3,i) 

X [l-f d'' (ri ,r2,t)+ d"^ (r2 , rg , i) 
+d^^'^nr3,ri,i)], (9) 

where the dislocation-dislocation correlation functions 
are defined in the usual way: 



Pi iri,t)pi^r2,t) 



(10) 



Substituting the cluster expansion ^ into the second- 
order evolution equation ([5]) leads to a closed evolution 
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equation for ^2^'**^, or equivalently, for d^^'"^. Its steady 
state solution was calculated by Vinogradov and Willis 
for infinite parallel screw dislocations.— Their motion 
was constrained to simple glide in one direction, so the 
difference between the system studied by them and the 
one considered in this paper is only the form of the gener- 
ated stress field of the particles Ti„d • They deduced an ex- 
act solution for the correlation function when there were 
only positive screw dislocations. However, in the case 
when both positive and negative screw dislocations were 
present, from the equation for the correlation functions 
they deduced a mathematical contradiction. Namely, the 
correlation functions had to tend to a constant nonzero 
value for large distances, which is obviously incorrect for 
real dislocation arrangements. They also showed that 
this contradiction was independent from the actual form 
of Tiiid , which means that for edge dislocations the equa- 
tions are not solvable either. We thus speculate that 
according to this result the cluster approximation is un- 
physical for the system of parallel edge dislocations. 

Another traditional method for closing the hier- 
archy at this order is the Kirkwood superposition 
approximation^'*, which was successfully used for in- 
stance in fiuid mechanics (for a brief overview see, e.g., 
Ref . [25h . In the context of dislocations, it was first intro- 
duced by Zaiser et al7 but after investigating the result- 
ing equations only a few basic conclusions were drawn. 

The Kirkwood superposition approximation expresses 
the three-particle density functions in terms of the pair 
densities: 

P3"'"''(n,r2,r3,i) 

'''{r,,r,, t)pl^ (r2 , r3 , Qp^ '^^ (^3 ,r,,t) 
pr(ri,t)p?^(r2,tVr(r3,0 
= pr(ri,t)pr(r2,iK(r3,i) 

X [1 + d^i'^^ (r-i, r2, + d'-^'^ {r2,rs, t)] 

x[l + d'''''{r3,ri,t)]. (11) 




si d 

and 




S2 d 



If we keep only the terms linear in the correlation func- 
tions d"^'"^ this expression gives the same result as the 
"cluster approximation" ([9]), i.e., their asymptotic be- 
havior for large distances is identical. But as we will see, 
the nonlinear terms play an important role in the forth- 
coming derivation. We would like to emphasize that the 
Kirkwood factorization is indeed an approximation. For 
example, it does not fulfill the simple normalization con- 
dition, that if the number of particles is finite (let's say 
N), then J^^ p3{., ., r, t)d'^r = {N - 2)p2. Thus it is ex- 
pected that the superposition approximation can't give 
precise results for small distances. 



In the following, we continue with substituting the ap- 
proximation (jlip into the second order evolution equa- 
tion ([5]). According to the previous results for the stud- 
ied homogeneous system [Eqs. ([5]) and the second 
term of ([8]) simplifies to 



= 2siS2 7^ [bp2'"'(n -r2,t)Tind(ri -r2)] . (12) 
or I 



For the third and the fourth terms of ^ one gets 



6p-'-(ri 



r2: 



P2'''"'i-'^3,t)pT'''^(r3 - ri + r2,t)Tind{ri - r2 ~ r^)d^r-i 



(13) 



S3=±l 



r2 



ri + r2,t)/92^''*'(-r3,i)riiid(T"i - ^2 - r3)d'^r3 



(14) 



where we also applied the T-mdi—f) = — 'nnd(^) relation. Finally, after introducing r .— ri — r2, for the evolution 
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equations of the second order one arrives at 



dt 



2siS2Tind(r) + (Sl - S2)Toxt + 



S3=±l 



= 0. (15) 



r 



The correlation functions, defined by Eq. (fTU]) . due to 
([5]) and (O simphfy to 



(16) 



By substituting into ITSl) one gets a closed set of 
equations for the two-particle correlation functions d^^^^^ . 
In the rest of this paper it is assumed, that: 

• The number of positive and negative signed dislo- 
cations are equal (p^ = p^). 

• The external shear stress is zero (Text = 0). 

As it is explained in details in the Appendix, it follows 
that d++{r,t) = d—{r,t) and d+-(r,t) = d-+{r,t) for 
every r and meaning there are only two independent 
correlation functions: and . Under these condi- 
tions the evolution equations simplify to (for details see 
the Appendix): 

dtd++ + 2V [6(1 + d++)(Ti„d + ri^ + + r^')] = 0, 

(17) 

dtd+- + 2V [6(1 + rf+-)(-ri„d - rj^ - + r^)] = 0, 

(18) 

where we have introduced the following terms having 
stress dimension: 



rtir^t) := p+ / 2dd{r' ,t)nUr - r')d\', 



(19) 



Tl^{r,t) ■.^p+ I 2dd(r',t)4(r-r',t)Tind(r-r')dV, 

(20) 

T^{r,t) :=p+ / 2ds{r',t)dd{r-r',t)n^d{r~r')d\' 



(21) 

with ds := id++ +d+-)/2 and d^ := (d++ -d+-)/2. 

Dislocations of opposite signs should move at each 
given point with equal velocities but in opposite direc- 
tions, which implies 



(22) 



for every r and t. This argument is identical to that in 
the appendix of Ref. where a similar stress term is 
omitted. 



The final evolution equations (|17p and (HH]) with 
are very similar to those obtained by Groma et al. for 
the one particle dislocation densities;^ To emphasize the 
analogy even more, let us introduce the following nota- 
tions: 

p'\r,t) ■.= p+ [2 + d++{r,t) + d+-{r,t)\ 

= 2p+[l + 4(r,i)], (23) 
K\r,t) ■.= p+ [d++{r,t)-d+-{r,t)] 

= 2p+dd{r,t). (24) 

(Here and previously the 'h' superscript refers to the ho- 
mogeneous system.) Although these quantities have den- 
sity dimensions, they are auxiliary quantities and do not 
carry the meaning of single particle densities. The evolu- 
tion equations for these newly introduced quantities can 
be written as 



dtp"" + 2V [6k" (rind + tI'^ + T^)] = 0, 
dtK^' + 2V [6p"(ri„d + t':, + T^)] = 



with 



T^Ar,t)= / K^r',t)Ti„d(r-r')dV, 



(25) 
(26) 



(27) 



T^^{r,t)^ U\r',t)ds{r-r',t)n^d{r-r')d\'. (28) 



These are the equations which govern the evolution of 
dislocation-dislocation correlation functions. Instead of 
their numerical solution, we will prove their correctness 
by an analogy with the phenomena of screening of an 
individual dislocation, which is discussed in the following 
section. 



III. DISLOCATION DENSITY FIELDS OF A 
SCREENED EXTERNAL DISLOCATION 

The stress field of a single dislocation decays as 1/r. 
This implies that if the distribution of the dislocations 
was completely random in a crystal, then the elastic en- 
ergy per unit volume would diverge logarithmically with 
the crystal sizci^ Since this cannot be observed in real 
systems, the only possible solution is, that dislocations 
arrange themselves in a correlated way which screens out 
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their long range effect. The phenomenon lends itself to 
the analogy with Debye screening in Coulomb systems. 

In order to address the problem, Groma et al^^ stud- 
ied the induced geometrically necessary dislocation den- 
sity around a single external edge dislocation in 2D, pro- 
posed an equation for the stress potential in equilibrium 
and gave analytic solution for the infinite plane, provided 
the total density was constant and much larger than the 
geometrically necessary dislocation density. While this 
is not the original question of the correlation function, 
within linear response theory, however, the result is ex- 
pected to be valid also to the correlation problem. It was 
indeed shown that the screened dislocation's stress field 
decays faster than that of the unscreened one's, in the 
direction perpendicular to the Burgers vector by a power 
law, and in all other directions exponentially. This re- 
sult was actually compared with correlation simulations, 
mostly because this was numerically simpler than the 
computation of the induced field by an external disloca- 
tion, and in the direction of the power decay not only the 
exponent but the entire shape predicted by theory was 
recovered in the simulation. However, the exponential 
decay along the x axis was not seen as in the theoretical 
solution, mostly because of the special conditions of dis- 
location motion on a torus as realized in the simulation. 
In any event, the equivalence of response to an external 
dislocation and the correlation in the absence of a fixed, 
external one, was taken as granted, precisely the problem 
addressed in the present paper. 

An important peculiar aspect of screening of disloca- 
tions has to be emphasized. First, we recall that in a 
thermal system like Coulomb plasma linear response the- 
ory provides enough ground to expect the similarity be- 
tween screening of an external charge and that in the 
correlation. On the other hand, the screening problem 
of dislocations already emerges at zero physical temper- 
ature. Now the role of the temperature in keeping oppo- 
sitely charged particles at distance is taken over by the 
constraint to slip axes, so here an effective temperature 
parameter arises that can be determined from compari- 
son with simulatioufSSiSi Now given the fact that we only 
have a temperature parameter (which directly appears 
in the effective thermodynamic potential) but not usual 
thermodynamics, it is by far not obvious that the two 
types of screening, namely, the one of an external dislo- 
cation, and the one appearing in the correlation function, 
should maintain the same type of equivalence as if Boltz- 
mannian thermodynamics were valid. Hence, in the case 
of dislocations it remains an open problem, what the re- 
lation is between screening by response to an external 
effect and screening in the correlations. 

In the previous section we constructed the equations 
of motion for correlations, now we do the same for the 
one-particle densities in the presence of an external dis- 
location. As a result of the stress field generated by the 
inserted object, the positions of the other dislocations 
change and a new relaxed state evolves. In it, for in- 
stance, one particle dislocation densities won't be con- 



stant any more. In this section, the time evolution of 
these functions is investigated. 

Due to the spatially varying extra force acting on the 
dislocations the system is not translation invariant any 
more. So, in the case of screening of an external dislo- 
cation, contrary to the case in the previous section, the 
system is spatially inhoniogeneous. The evolution of the 
one-particle densities is described by the first member of 
the BBGKY hierarchy HI): 

dt 

S2=±l„2 

(29) 



+^b^(n,t)f;^Un)] = o, 

where we introduced the notation F^^^.{r) := s6bTi„d('") 
with b := \b\ for the stress field of the extra dislocation. 

According to the deduction in Ref. ^, Eq. ((29|) can be 
cast into the following form 

dtp + V [bK{Tind + Tsc -Tf+ Tfc)] = 0, (30) 
dtK + V [bp{Tind + Tsc -Tf+ Tfc)] ^ 0, (31) 

where we have introduced the 

p{r,t):^p+{r,t)+p-{r,t) (32) 

and 

K{r,t):=p+{r,t)~p-{r,t), (33) 

total and geometrically necessary dislocation densities, 
respectively, and 



T,e(r,t) / K{r',t)Tind{r-r')d\', 



(34) 



n{r,t) 



Tf{r,t) 



K{r', t)d{r, r', i)ri„d(r- - r')d^r', (35) 



p{r',t)da{r,r',t)nnd{r-r')d\' 



(36) 



with d := {d++ + d — + d+^ + d"+)/4 and da := 
(d^ — d ^)/2. We note that in this inhomogeneous 
case the correlation functions do not depend on the dif- 
ference of their arguments, meaning one has to remain at 
the general definition (fTO| . Here and in the rest of this 
paper (~) indicates that the correlation function has two 
spatial variables. During the derivation of Eqs. ([501 EJ) 
no approximations have been made. 

If the correlation functions were short-range, it could 
be approximated that they depend only on the relative 
coordinate ri — 

Pl'^'\r,,r2,t) 

= (n , i)pf (r2 ,t)[l + (n - r2 , t)] . (37) 
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Here the correlation function d'^^'"^ can be taken 
from homogeneous systems The assumed shortness of 
dislocation-dislocation correlations was proved on dis- 
crete dislocation simulation results earlier?^ It was found 
that the correlation functions decay to zero exponentially 
within a few average dislocation spacings. 

Because of the short-range nature of dislocation- 
dislocation correlations, the k and p functions can be 
approximated by their Taylor expansions in the integrals 
of Eqs. ([55)1 and ([55)1 . After keeping only the first non- 
vanishing terms Groma et al. arrived at 

^bir,t)^-—^ -Dd-^d,Kir,t) (38) 



2^(1 -z.) "p(r,t) 



and 



rfir, t) = —^i^CrffevMM), (39) 

where Dd and Cd are positive constants.* The term Tf, 
is a gradient like contribution to the stress, and is called 
"back-stress" , while r/ can be interpreted as a local flow- 
stress. The physical correctness of these approximations 
was proved by discrete dislocation simulations 

The description of Debye screening of an external dis- 
locations was based on a thermodynamic potential-like 
functional, from which the evolution equations ((30l [3T|) 
for K and p can be derived . In it Tf, is approximated 
as (|55)) and t/ is omitted. In the case of p = const, 
the analytical solution of the static case was given and 
compared with numerical simulations ^ 

In what follows we shall see that, to establish the con- 
nection between the evolution of correlation functions 
and the screening field of a single external dislocations, 
we will not need the approximations ([571) . (ES), and ((55)) 
for the correlation functions, t;,, and Tf, respectively, 
rather can keep the more general definitions (|35p and 
(ESI). 



IV. CONNECTION BETWEEN THE 
SCREENED DISLOCATION DENSITIES AND 
THE CORRELATION FUNCTIONS 

In the previous two sections we have derived evolu- 
tion equations for correlation functions in a homogeneous 
system, within the Kirkwood approximation, as well as 
for induced dislocation densities around a fixed, exter- 
nal, dislocation. A main message of this paper is that 
the deduced equations for the two different cases, namely 
(I23I24P and (|30l3ip . are very similar. More is true, how- 
ever, they become identical, if special equalities hold for 
various d-s and time is appropriately rescaled. In detail, 
these conditions arc the following: 

• If the numbers of + and — dislocations are the 
same, then in the absence of external stresses 
d+-{r,t) = d~+{r,t) and d++{r,t) = d--(r,t). 



So, if the general correlation functions are approx- 
imated by the ones taken from homogeneous sys- 
tems, i.e.. 



(40) 



then da{ri,r2,t) = and d(ri,r2,t) = ds{ri - 
r2,t). In this case the "fiow-stress" term Tf dis- 
appears from pOI3ip . while it was never present 
in Eqs. (gSl [M]). Moreover, the Tb in Eqs. (|30I31|) 
becomes equal to the of (|23l24p . 

The main theoretical result of this paper immedi- 
ately follows. Namely, if one takes a solution of 
^+(r,t) and d^ (^i^); '^ext sub- 
then 



Eqs. im^2^ 
stitutes these 



for d^ 
into 



30\31\) . then the latter equa- 
tions will have a solution identical with the one of 

fug). 



Physically, the equality d{ri,r2,t) = ds{ri — r2,t) 
corresponds to the approximation when the effect 
of the stress by the external dislocation in Eqs. 
(|30l3ip are neglected in the correlation functions. 
This condition amounts to taking the external ef- 
fect as a first order perturbation, in perfect anal- 
ogy of linear response theory valid for conventional 
thermal systems. 



• The divergence-like terms of ([25)) and (|M)) are mul- 
tiplied by two, which only affects the time scale 
of the process. The appearance of this factor can 
be well understood if we consider two dislocations 
moving in each other's stress fields. They relax into 
the same configuration either we keep the position 
of the first dislocation fixed or not. The only differ- 
ence is, that the relaxation in the first case (which 
corresponds to screening of an individual disloca- 
tion) lasts twice as long as in the case when both 
dislocations can glide freely (case of the evolution 
of correlation functions). 

To sum up, since a.t t — p{r, 0) = /o''(r, 0) = /?+ + p^ 
and K{r,0) = K'^{r,0) — 0, it follows that 



p{r, 2t) = p\v, t) =p+[2 + d++{r, t) + d+-{r, t)] 



(41) 



and 



K(r,2i) = K''(r,t) ^ p+ [d++{r,t) - d+-{r,t)] (42) 

for every r and t. We would like to stress, that during 
the derivation of this statement we made approximations 
for both {p^,K^) and for {p,k) functions. For the first 
set the Kirkwood superposition approximation [Eq. pip ], 
and for the second linear response was assumed, mean- 
ing general correlation functions were approximated by 
ones taken from infinite homogeneous systems [Eq. (j^O)) ]. 
If the comparison of the numerically obtained functions 
confirms the validity of (|4T)) and (|42)) . then it would also 
confirm the applicability of the Kirkwood superposition 
approximation. In the following section this comparison 
is discussed. 
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SIMULATION RESULTS 



following function was considered: 



To measure p'', p, k^, and k numerically, we first per- 
formed more than 13000 different relaxations. In each of 
them at the beginning 128 dislocations (64 with positive 
and 64 with negative sign) were distributed randomly 
in an L X L square-like domain with periodic boundary 
conditions, and then they relaxed to a steady state. The 
and correlation functions and then the p'^ and 
quantities were determined from these configurations 
by counting the relative coordinates of the dislocations. 
They are plotted in the first column of Fig. [TJ 

Afterwards, a new dislocation with positive sign was 
placed at the center of the simulation area in each con- 
figuration, and the systems were let to relax again while 
the extra dislocation was fixed. These simulations corre- 
spond to the phenomena of Debye screening of an exter- 
nal dislocation. The p+ and densities were obtained 
by counting the number of positive and negative disloca- 
tions falling into each cell of a 256 x 256 mesh. Then p 
and K functions defined by (|32p and ((33)) were calculated. 
They are also plotted in Fig. [T] 



\p{r)-dp''{r)\d^r. 



(43) 




-3-2-10123 



-3-2-10123 



FIG. 1: The p , p, k , and k functions in the relaxed 
state obtained from numerical simulations. The distances 
are measured in pj^?'^ average dislocation spacings, where 
ptot := P^ + P~ is the total dislocation density. 



The two functions, according to our expectations, are 
quite similar. The difference is that in the case of k and 
p the amplitude of the noise is much higher. The simple 
reason for this is that a relaxed system of N dislocations 
gives N{N — I) data (the relative coordinates) during 
the calculation of the correlation functions but only N 
data (the positions) when calculating the densities. In 
order to make a comparison between these functions the 




The integral of a noise is nearly constant so if p and p^ are 
equal except for the noise then this quantity must have 
a minimum at d = 1. As it was mentioned in Sec.[Tll the 
Kirkwood superposition approximation is not expected 
to be correct for distances smaller than an average dislo- 
cation spacing, and so, Eqs. (|¥T|) and should apply 
only for r = |r| > p^?-^ (ptot := P+ + P' and p^?'^ is 
the average dislocation spacing). Accordingly, it is mean- 
ingful to define the domain of the integration A so that 
a circle of radius 0.5p^^'^ around the origin is left out 
from it. The numerically obtained function Ip is plotted 
in Fig. [21 As it can be seen, Ip indeed has a minimum 
around d = 1. This confirms that the p^ and p fimctions 
can be considered equal in first approximation. The rep- 
etition of this calculation for and n leads to similar 
results. 



17 I 
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FIG. 2: The Ip function calculated numerically at different 
values (data points). A second order polynomial was fitted 
around the minimum of the points (solid line). 

In Fig. [3] p'* and p functions can be seen at different 
times obtained from the numerical simulations. Accord- 
ing to (|4T|l . p''(r, t) was compared with p(r, 2t), because 
of the factor of 2 appearing in the evolution equations of 
the correlation functions in Eqs. ([^S)) and The same 
can be seen for the k functions in Fig. |4l To summarize, 
it can be stated that (|4l]l and (I42|) are, at least in first 
order, fulfilled. 



VI. CONCLUSIONS 

In this paper we derived evolution equations for the 
correlation functions of long parallel edge dislocations. 
We started our analysis at the BBGKY hierarchy of dis- 
location many-body densities, which was deduced earlier 
from the equations of motion of individual dislocations)^ 
As a closure approximation for the hierarchy, the Kirk- 
wood superposition approximation was used. Our in- 
vestigations showed that the evolution of the correla- 
tion functions and the dislocation density fields evolv- 
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FIG. 3: The evolution of (first row) and p (second row). The time t is measured in 
B is the dislocation mobility. According to Eq. (|41|l p'^{r,t) and p{r,2t) are compared. 



b-'fj. 



2-k{1-v)BL' 



^ dimensionless units, where 
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FIG. 4: The evolution of (first row) and k (second row). The time t is measured in 2-w(i-v)bl'^ dimensionless units, where 
B is the dislocation mobility. According to Eq. (I42|l K*^{r,t) and K{r,2t) are compared. 



ing around a screened dislocation are closely related (ac- 
cording to Eqs. (|4T|) and (|42)) only the time scale of the 
processes is different). In Sec.|V]this statement was con- 
firmed by numerical simulation. 

The Kirkwood superposition approximation pT|) gives 
the three-particle dislocation density in terms of the pair 
correlations. Contrary to the cluster approximation ([9]), 
it contains quadratic and cubic terms of the correlation 
function beside the linear one. As it was mentioned in 
Sec. In] these higher-order terms play an important role. 



According to our analysis presented in this paper, the 
omission of the cubic term leads to the disappearance of 
the back-stress term in the evolution equations P5|) 
and (|26p . This corresponds to the mean- field approxima- 
tion where dislocation correlations are neglected result- 
ing that the stress acting on dislocations is simply the 
self-consistent field t^^^ 

If one neglects the cubic and quadratic terms (which 
is the cluster approximation), not only the back-stress 
T^ , but also the self-consistent field r^'j. disappears, which 
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is obviously unphysical. We mentioned in Sec. |TT] that 
the cluster approximation ^ used by Vinogradov and 
Willis leads to a mathematical contradiction in the case 
of edge dislocations of opposite signs. ^^'^'^ This result is 
in complete agreement with our analysis. 

The Kirkwood superposition approximation leads to 
the evolution equations and P5|) . Their correctness 
was proved implicitly by the facts that (i) analogy was 
found with the one-particle density evolution equations 
and (ii) numerical work confirmed the similitude between 
the evolution of the one- and two-particle distributions. 
Another way of testing the deduced integro-diffcrcntial 
evolution equations could be their numerical solution. 
This is out of the scope of the present paper, and re- 
quires further numerical work. 
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APPENDIX 

In this section we derive evolution equations for the 
correlation functions [which were defined by Eq. PTI)) ] 
starting from the evolution equation lfT5|) of jOj^'"^. First, 
we assume that the densities of positive and negative 
dislocations are equal {p'^ = p~). It can be noted that 
in this case if p^+(r,i) = P2^{r,t) for any t, then the 
second term of ([T5|) is equal for both and P2 - This 
means, that the evolutions of and P2~ are identical, 
hence P2^{r,t) = p2~{r,t) and therefore d++(r,t) = 
d for every t. For this reason one only has to 

deal with d^^ , d"' , and d ^ . 

Let us introduce 



and 



dp := {d+- +d-+)/2 



da id+-^d-+)/2, 



(A.l) 



(A.2) 



namely the symmetric and the antisymmetric part of 
d"' . Furthermore, we define dg and dd as the half of 
the sum and the difference of d^'^ and dp, respectively: 



and 



id+++dp)/2 



dd := (d++ - dp)/2. 



(A.3) 
(A.4) 



It is also useful to introduce the following quantities hav- 
ing stress dimensions: 



T,',{r,t) :=p+ / 2d<i(r',t)rind(r-r')dV', 



(A.5) 



Tl^{r,t) ■.= p+ / 2dd(r',04(r-r',i)Tind(r-r')dV, 



(A.6) 



T>f{r,t) ■.= -\p+ I da{r',t)dair-r',t) 



R2 



(A.7) 



2dsir',t)dd{r-r',t) 



+ -da{r',t)da{r~r',t) 



T,^d{r - r')d\\ (A., 
rp^(r,i) :=p+ / \ da{r\t) [l + d++ {r - r\t)\ 



- [1 + d++(r', t)] da{r - r',<)}Ti„d(r - r')d^r' . 

(A.9) 

With these new functions, from the evolution equations 
(fT5|) for the correlation functions one obtains: 

dtd++ + 2V[b{l + d++) 

X (Tind + ri^ + - t'} + t':)\= 0, (A. 10) 

dtd+^ +2V[b{\+d+-) 

X (-Tind + Text - t!}^ ~ + + + T^)] = 0, 

(A.ll) 

dtd-+ + 2\/[b{l+d-+) 

X (-Tind - Text - T,", - T,^ + + - T^)] = 0. 

(A.12) 

This result can be checked by simple substitution. 

In the absence of external stress {Text — 0) the gen- 
eral equations (jA.lOl - IA.12P can be simplified because 
for symmetry reasons d^ (''j^) = d ^{r,t) for every t. 
In other words, at zero external stress if one alters the 
sign of all dislocations it does not affect the behavior of 
the system. This means that da = and so rj? = = 0. 
In this case the resulting evolution equations are: 

dtd++ + 2V [b(l + d++)(Ti„d + + + T^)] - 0, 

(A.13) 

dtd+- + 2V [6(1 + d+-)(-ri„d - ri'c - + = 0- 

(A. 14) 
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